Fourth-order Correlations of Conserved Charges in the QCD 

Thermodynamics 

Wei-jie Fu^Q and Yue-liang Wu 1 '^ 

1 Kavli Institute for Theoretical Physics China (KITPC), 
Key Laboratory of Frontiers in Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Science, Beijing 100190, China 
(Dated: October 6, 2010) 

Abstract 

The fourth-order correlations of conserved charges, such as the baryon number, electric charge, 
and strangeness, are studied at finite temperature and nonzero baryon chemical potential in an 
effective model. It is found that the fourth-order correlations change rapidly and have three extrema 
during the chiral crossover with the increase of the baryon chemical potential. The absolute values 
of the extrema approach infinity when the thermodynamical system moves toward the QCD critical 
point and the fourth-order correlations are divergent at the critical point. The contour plots of the 
fourth-order correlations in the plane of the temperature and baryon chemical potential are given. 
It is noticed that all the fourth-order correlations of conserved charges, except for X13 an d X13 \ 
are excellent probes to explore the QCD critical point in heavy ion collision experiments. 
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It is believed that the deconfined quark gluon plasma (QGP) is produced in ultrarela- 
tivistic heavy ion collisions jl-8|. Therefore, QCD thermodynamics, such as the equation of 
state of the QGP, chiral and deconfinement phase transitions, QCD phase diagram and so 
on, has been a subject of intensive investigation in recent years. One distinct characteristic 
of the QCD thermodynamics is that there is a critical point in the QCD phase diagram 
in the plane of temperature and baryon chemical potential, which separates the first-order 
phase transition at high baryon chemical potential from the continuous crossover at high 
temperature This characteristic is confirmed in various field theory models 10-18]. Al- 
though there is no definite evidence that the QCD critical point also exists in the lattice 
QCD calculations due to the sign problem at finite chemical potential, some lattice groups 
find that the QCD critical point maybe exist in the phase diagram, based on extrapolating 
results at small Ub/T (ratio of the baryon chemical potential and the temperature) to those 
at large /j^b/T [19|-|21|. In the meantime, experiments with the goal to search for the QCD 
critical point are planned and underway at the Relativistic Heavy Ion Collider (RHIC) at 



the Brookhaven Nationa 
CERN in Geneva 



22 



25]. 



Laboratory (BNL) and at the Super Proton Synchrotron (SPS) at 



Then, searching for and locating the QCD critical point becomes a crucial and vital task. 
It has been proposed that the QCD critical point can be found through the non-monotonic 
behavior of fluctuations and correlations of various particle multiplicities as functions of 
varying control parameters 26M36| . Particularly, fluctuations and correlations of conserved 
charges, such as the baryon number, electric charge, and strangeness, deserve more atten- 
tions. On the one hand, the fluctuations and correlations of conserved charges are sensitive 
to the structure of the thermal strongly interacting matter and behave differently between 



the hadronic and QGP phases 
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27, 3l|, |37j. On the other hand, since the conserved 



charges are conserved through the evolution of the fire ball, the fluctuations and correlations 
of conserved charges can be measured in heavy ion collision experiments. 
In our previous work 



, we have calculated the fluctuations of conserved charges up to 
the fourth-order and the correlations to the third-order at finite temperature and nonzero 
baryon chemical potential in the 2+1 flavor Polyakov-Nambu-Jona-Lasinio (PNJL) model. 
We found an interesting result: among all the fluctuations and correlations discussed in our 



previous work 



the numerical calculations indicate that Xn > X21 \ X21 1 an d Xnf (those 



notations will be introduced in the following) are the most valuable probes for exploring the 



2 



QCD critical point. All these quantities are the third-order correlations of conserved charges. 
Therefore, it is natural to go beyond our previous work to study the fourth-order correlations 
of conserved charges near the QCD critical point. Furthermore, using higher-order moments 
of particle multiplicity distributions to search for the QCD critical point is possible in the 
Beam Energy Scan at RHIC {24]. In this work, we will calculate the fourth-order correlations 
of conserved charges at finite temperature and nonzero baryon chemical potential in the 
PNJL model. Most attentions will be paid on the studies of the non-monotonic behavior of 
the fourth-order correlations near the QCD critical point. 

We begin with the definition of the correlations of conserved charges as follows 

BQS _ d ^+ k (P/T*) 

Xi i k " d(f, B /Tyd^ Q /Tyd^ s /Tr 1 1 

where P and T is the pressure and temperature of a thermodynamical system; ^b,q,s are 
the chemical potentials for baryon number, electric charge, and strangeness, respectively. 
These conserved charge chemical potentials are related with the quark chemical potentials 
through the following relations, 

12 11,11 

fa = + g/XQ, /id = g/ifl - g/XQ, and fx 8 = -hb - g/xc? - Ms- (2) 

where fi Ut d, s are the chemical potentials for u, d, and s quarks, respectively. Denoting the 
ensemble average of conserved charge number Nx {X = B, Q, S) with (Nx), we can obtain 
the fourth-order correlations as follow 

X? 3 y = ^((6N X 5N Y 3 ) - 3(6N x 6N Y )(6N Y 2 )y (3) 
x! Y = ^r 3 ((SN x '6N Y 2 )-(5N x 2 )(6N Y 2 } 

-2(5N X 5N Y ) 2 ^, (4) 
X21V = ^[(SN X 2 5N Y 5N Z ) -2(5N X 5N Y )(5N X 5N Z ) 

-(5N x 2 )(5N Y 6N z }y (5) 

where 8N X = Nx — (Nx) and V is the volume of the system. 

We adopt the 2+1 flavor Polyakov-loop improved NJL model to study the fourth-order 
correlations of conserved charges near the QCD critical point. The validity of this effective 
model is expected, since the critical behavior of the QCD phase transition is governed by 
the universality class of the chiral symmetry, which is kept in this model. Furthermore, 



compared with the conventional Nambu-Jona-Lasinio model, the PNJL model not only has 
the chiral symmetry and its dynamical breakin g m echanism, but also includes the effect 
of color confinement through the Polyakov loop 18|, |39H46| . Furthermore, the fluctuations 



flavor PNJL model at finite 



36[ are well consistent with 



and correlations of conserved charges calculated in the 2+1 
temperature but with vanishing chemical potentials in Ref. 
those obtained in lattice calculations [47], which shows that the 2+1 flavor PNJL model is 
well applicable to study the cumulants of conserved charge multiplicity distributions. 
The Lagrangian density for the 2+1 flavor PNJL model is given as 181 ] 
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£ PN jl = ^(hii D>1 + 7oA - m )tp + G ^2 [(V^) 2 + {$i>l$T a ip)' 



a=0 



—K 



det/Ml + 7 bM + det/$(l - - $* , T) 



where ip = (ip u , ipd, r 'Ps) T is the three-flavor quark field, and 



iA" with A" = 5qA 



A u 



gA ° a Y 



-iA, 



4- 



(6) 



(7) 



where A a 's are the Gell-Mann matrices in color space and g is the gauge coupling strength. 
rho = diag(mQ, ttiq, mg) is the three-flavor current quark mass matrix. Throughout this 



work, we take mj( 



mo, while keep m s Q being larger than m , which breaks the 



577(3)/ symmetry. The matrix p, = diag(/x u , /x^, /x s ) denotes the quark chemical potentials 
which are related with the conserved charge chemical potentials through relations in Eq.(j2]). 

In the above PNJL Lagrangian, £/($,$*, T) is the Polyakov-loop effective potential, 
which is expressed in terms of the traced Polyakov-loop $ = (Tr c L)/N c and its conjugate 
$* = (Tr c tf)/N c with the Polyakov-loop L being a matrix in color space given explicitly by 



L(x) = V exp i / dr A±(x,t) = exp[i/M 4 



(8) 



with (3 = 1/T being the inverse of temperature and A^ = iA°. 

In our work, we use the Polyakov-loop effective potential which is a polynomial in $ and 



421 ]. given by 



with 



T 4 2 6 v 1 4 v ; 



6 2 (T) = a + ai(^) + a 2 (^) + a 3 (^) . 



(9) 
(10) 



4 



The parameters in the effective potential are fitted to reproduce the thermodynamical be- 
havior of the pure-gauge QCD obtained from the lattice simulations, and their values are 
a$ = 6.75, ai = —1.95, a<i = 2.625, a 3 = —7.44, 6 3 = 0.75, and 64 = 7.5. The parame- 
ter T is the critical temperature for the deconfinement phase transition to take place in 
pure-gauge QCD and T Q is chosen to be 270 MeV according to the lattice calculations. As 
for the five parameters in the quark sector of the model, we adopt their values which are 
ml = 5.5 MeV, m s = 140.7 MeV, GA 2 = 1.835, KA 5 = 12.36, and A = 602.3 MeV, which 



are fixed by fitting t 
and f w = 92.4 MeV 



re observables m n = 135.0 MeV, m K = 497.7 MeV, m v > = 957.8 MeV, 
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In the mean field approximation, the thermodynamical potential density (Q = — P) for 
the 2+1 flavor quark system is given by 

f=u,d,s 

+ Z l n [1 + Z<&* e -i.4-Pf)/T + 2$ e -2(E*-p f yT + e -3(E/-^)/Tj 

3 

+ I. l n [l _|_ 3$ e -(^+M/)/T + 3$* e -2(^+M/)/T + e -3(J^+M/)/Tj | 

3 J 

+2G(0 U 2 + <p d 2 + S 2 ) - AKct> u <f> d ct> s + $*, T), (11) 

where fa's (i = u, d, s) are the quark chiral condensates, and the energy-momentum 
dispersion relation is E l p = ^/p 2 + M 2 , with the constituent mass being 

M { = ml - 4G<f) t + 2K<f )j </> k . (12) 

Minimizing the thermodynamical potential in Eq. ( II ip with respective to <ft u , (pa, <Ps, and 
$*, we obtain a set of equations for the minimal conditions, which can be solved as functions 
of temperature T and three conserved charge chemical potentials [Iq, and ^5. 

We use the method of Taylor expansion to compute the fourth-order correlations of 
conserved charges in the PNJL model. In Fig. [1] we show the fourth-order correlations 
between two conserved charges versus the baryon chemical potential at several values of 
temperature calculated in the PNJL model. We find that the QCD critical point is located at 
about T c = 160 MeV and ^b c — 819 MeV (/iq = [As — 0) with input parameters given above. 
From Fig. HJ one can clearly see that the magnitudes of all the fourth-order correlations in 
this figure grow rapidly and oscillate drastically when the QCD phase transition occurs. 
Comparing different curves corresponding to different temperatures, we can recognize that 
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FIG. 1: (color online). Fourth-order correlations X31 (top- left), X22 (top-middle), X13 (top- 
right), X31 (center-left), X22 (center-middle), X13 (center-right), X31 (bottom-left), X22 (bottom- 
middle), and x% (bottom-right) as functions of the baryon chemical potential \ib (hq = /As = 0) 
with several values of temperature in the PNJL model. 

the oscillating amplitudes of the correlations increase quickly when moving toward the QCD 
critical point, i.e., the temperature approaches to T c = 160 MeV. All the fourth-order 
correlations between two conserved charges diverge at the QCD critical point. When the 
temperature is below T c , the chiral phase transition is first-order and the correlations are 
discontinuous during the phase transition as shown by the dashed and dotted lines in Fig. [TJ 
When the temperature is above T c , the chiral phase transition is a continuous crossover due 
to finite quark current mass. Correspondingly, the correlations of conserved charges are also 
continuous during the QCD phase transition. As for the fourth-order correlations between 
two conserved charges, one notices that there are two maxima and one minimum in the 
curves of X31 7 X22 7 X13 7 X22 7 an d X22 > while there are two minima and one maximum in 
the curves of other fourth-order correlations between two conserved charges. Furthermore, 
it is seen that all the fourth-order correlations except for xW an d X^z vanish rapidly once 
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the thermodynamical system deviates from the QCD phase transition. In Fig. [I] we observe 
that X13 an d X?3 still have finite values at large baryon chemical potential, this is because 
the strange quark have relatively large current mass and its constituent mass can not be 
neglected even in the chiral symmetric phase. 




FIG. 2: (color online). Contour plots of the fourth-order correlations X31 (top-left), X22 (top- 
middle), X13 (top-right), xif (center-left), X22 (center- middle) , X13 (center-right), X31 (bottom- 
left), X22 (bottom-middle), and X\s (bottom-right) versus temperature T and baryon chemical 
potential [Ib (fJ>Q = us = 0) in the PNJL model. 

Fig. [2] shows the contour plots of the fourth-order correlations between two conserved 
charges as functions of the temperature and baryon chemical potential calculated in the 
PNJL model. We observe that the fourth-order correlations between two conserved charges 
are vanishing in the chiral symmetry broken phase, i.e., in the bottom- left region of every 
little plot. It is noticed that the chiral phase transition line is distinct in these plots and one 
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FIG. 3: (color online). Fourth-order correlations among three conserved charges Xm^ (top), Xm^ 
(middle), and X211 (bottom) as functions of the baryon chemical potential //_b (a*q = /As = 0) 
with several values of temperature in the PNJL model. 



can easily recognize the region near around the QCD critical point, where the contour lines 
are dense. Comparing all these fourth-order correlations between two conserved charges, 
we notice that xfi Q , X22 \ X13 \ X31, X22 \ Xsf , and X22 are superior to Xrf and Xrf to 
be used to search for the critical point in heavy ion collision experiments, since the former 
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seven correlations have finite values only when the thermodynamical system is near around 
the QCD critical point. 
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FIG. 4: (color online). Contour plots of the fourth-order correlations among three conserved 
charges X\2i S (top), Xm (middle), and (bottom) versus temperature T and baryon chemical 

potential [ib (l^Q = fJ-s = 0) in the PNJL model. 

In Fig. [3] we present the fourth-order correlations among three conserved charges as 
functions of the baryon chemical potential with several values of temperature calculated in 
the PNJL model. It is seen that the fourth-order correlations among three conserved charges 
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are divergent at the QCD critical point same as those between two conserved charges. When 
the temperature is below T c = 160 MeV, those correlations are discontinuous at the first- 
order chiral phase transition; while when the temperature is above T c , they are continuous 
functions of the baryon chemical potential and change rapidly at the chiral crossover. It is 
observed that there are two minima and one maximum in the curves of Xi2i S an d X^n > 
and two maxima and one minimum in the curve of Xiu ■ O ne can a l so clearly notice that 
the fourth-order correlations among three conserved charges approach zero rapidly once the 
thermodynamical system deviates from the chiral phase transition. 

In Fig.H]we show the contour plots of the fourth-order correlations among three conserved 
charges as functions of the temperature and baryon chemical potential in the PNJL model. 
One can clearly recognize the QCD critical point in these three plots, which demonstrates 
that the fourth-order correlations among three conserved charges are quite sensitive to the 
singular structure related to the critical point. Therefore, the fourth-order correlations 
among three conserved charges are excellent probes to explore the QCD critical point in 
heavy ion collision experiments. 

In summary, we have studied the fourth-order correlations of conserved charges, such 
as the baryon number, electric charge, and strangeness, at finite temperature and nonzero 
baryon chemical potential in the 2+1 favor Polyakov-loop improved Nambu-Jona-Lasinio 
model. It is found that the fourth-order correlations of conserved charges are divergent at 
the QCD critical point. When the temperature is below the critical temperature of the QCD 
critical point, the fourth-order correlations are discontinuous at the first-order chiral phase 
transition; while when the temperature is above the critical temperature, the fourth-order 
correlations of conserved charges are continuous with the change of the baryon chemical 
potential and oscillate rapidly at the chiral crossover. In the curves of the fourth-order 
correlations as functions of the baryon chemical potential, there are three extrema during 
the chiral crossover. Furthermore, we have given the contour plots of the fourth-order 
correlations of conserved charges in the plane of temperature and baryon chemical potential. 
The QCD critical point can be easily recognized in these plots. Comparing the fourth-order 
correlations, we notice that all the fourth-order correlations of conserved charges, except for 
Xii and X13, are excellent probes to explore the QCD critical point, since these fourth- 
order correlations approach zero rapidly once the thermodynamical system deviates from 
the QCD critical point. 
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